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ABSTRACT 

The Kuiper Belt is a remnant from the early solar system and its size distri- 
bution contains many important constraints that can be used to test models of 
planet formation and coUisional evolution. We show, by comparing observations 
with theoretical models, that the observed Kuiper Belt size distribution is well 
matched by coagulation models, which start from an initial planetesimal popu- 
lation with radii of about 1 km, and subsequent coUisional evolution. We find 
that the observed size distribution above i? ~ 30 km is primordial, i.e., it has 
not been modified by coUisional evolution over the age of the solar system, and 
that the size distribution below ~ 30 km has been modified by collisions and 
that its slope is well matched by coUisional evolution models that use published 
strength laws. We investigate in detail the resulting size distribution of bodies 
ranging from 0.01 km to 30 km and find that its slope changes several times as a 
function of radius before approaching the expected value for an equilibrium coUi- 
sional cascade of material strength dominated bodies for i? < 0.1 km. Compared 
to a single power law size distribution that would span the whole range from 
0.01 km to 30 km, we find in general a strong deficit of bodies around ~ 10 km 
and a strong excess of bodies around 2 km in radius. This deficit and excess of 
bodies are caused by the planetesimal size distribution left over from the run- 
away growth phase, which left most of the initial mass in small planetesimals, 
while only a small fraction of the total mass is converted into large protoplan- 
ets. This excess mass in small planetesimals leaves a permanent signature in the 
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size distribution of small bodies that is not erased after 4.5 Gyrs of collisional 
evolution. Observations of the small KBO size distribution can therefore test if 
large KBOs grew as a result of runaway growth and constrain the initial plan- 
etesimal sizes. We find that results from recent KBO occultation surveys and the 
observed KBO size distribution can be best matched by an initial planetesimal 
population that contained about equal mass per logarithmic mass bin in bodies 
ranging from 0.4 km to 4 km in radius. We further find that we cannot match the 
observed KBO size distribution if most of the planetesimal mass was contained 
in bodies that were 10 km in radius or larger, simply because their resulting size 
distribution cannot be sufficiently depleted over 4.5 Gyrs to match observations. 

Subject headings: circumstellar matter — Kuiper belt: general — planetary 
systems — protoplanetary disks — planets and satellites: formation — minor 
planets, asteroids: general 



1. INTRODUCTION 



The Kuiper Belt consists of a disk of icy objects located just beyond the orbit of Neptune. 
In the Kuiper Belt planet formation never proceeded all the way to completion, which makes 
it an ideal laboratory for testing planet formation theories. 

The Kuiper belt size distribution contains many important clues concerning the forma- 
tion of Kuiper belt objects (KBOs), their effective strength, and their collisio nal evolution 
flDohnanvil[l969[ bavis fc Farinellal[l997i iKenvon fc Luul[l999[ |Pan fc Sar]l2005h . The cumu- 
lative size distribution of KBOs larger than i? > 30 km (i.e., objects with R-band magnitudes 
brighter than about 25) is well described by a single power-law given by 

N{> R) oc B'-'i (1) 

where A^(> R) is the number of objects with radii greater than i?, and q is the power-law 
index. Kuiper belt surveys find that the size distr ibution for KBOs with radii greater than 
about 30 km follows this power-law with q ~ 4 (e.g. lTrujillo et al.ll2001uBernstein et al.ll2004j : 
Fuentes fc Holmarj|2008l : iFraser et al.ll2008l ). which implies roughly equal mass per logarith- 
mic mass interval. This size distribution is a relic of the accretion history in the Kuiper belt 
and therefore provides valuable insights into the formation of large KBOs (R > 30 km) (e.g. 
SternlllQod : lOavis fc Farinellalll997t iKenvon fc BromlevI [20o3 ISchhchting fc Saril l201ll ) . It 
has been shown by several works that the la rge KBO size d i stribution can be well matched 



by numerical coagulation simulations (e.g. iKenyon fc Luul Il999l : ISchlichting fc Saril 12011 
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Kenyon fc Bromleyl |2012| ). For example, ISchlichting &: Saril (j201l[ ) find that the size dis- 



tribution of larger KBOs is well matched by planet formation models of runaway growth. 
During runaway growth only a small fraction of the total mass is converted into large proto- 
planets, while most of the initial mass remains in small planetesimals. The size distribution 
of the large protoplanets in the runaway tail follows a power law size distribution with dif- 
ferential power law index g ~ 4, implying roughly equal mass per logarithmic mass bin (see 
Figure [1]) . 




Radius [km] 

Fig. 1. — The size distribution at the end of runaway growth before the onset of coUisional 
erosion is given by the thick blue line. Note that during runaway growth, most of the initial 
mass remains in small planetesimals, while a small fraction of the total mass is converted 
into large protoplanets/KBOs. This specific example corresponds to an initial planetesimal 
population of bodies that were all 1 km in radius. As shown in Figures [2] and [3l the current 
large KBO size distribution is well matched by the resulting size distribution from runaway 
growth. For comparison, a power-law size distribution with differential power-law index 
q = 4 is given by the thin black line. 



Figures |2] and [3] show a dire ct comparison between t he results of runaway growth 
from the coagulation model from ISchlichting &: Saril (l201ll ) and the observed size distri- 
bution of dynamically cold and hot KBOs, respectively. The observed KBO size distri- 
bution was derived by iFuentes et al.l (120101 ) by combining results from KBO surveys by 
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cold refers here to objects with inclinations less than 5°, whereas dynamically hot corre- 
sponds to t hose with i > 5°. T h e agre ement of the observations with the simple coagulation 
model from ISchlichting fc Saril ( 1201 ll ) is good. Figures [2] and [3] show that the cold and hot 
populations can both be fit by the same size distribution, with the notable difference that 
the largest bodies in each population grew to different typical sizes. 

Provided that the fall off at large KBO sizes is not due to some selection effect, this 
suggests that KBOs grew to typical radii of about 100 km in the cold population and to 
typical radii of about 300 km in the hot population. If the growth in the hot and cold 
populations was terminated simultaneously, presumably by the excitation of the velocity 
dispersion of the growing KBOs and the smaller planetesimals, then this suggests that the 
hot population may have formed closer to the Sun than the cold population, because the 
shorter orbital periods and likely higher mass surface densities ensure faster growth at smaller 
semi-major axis. 

Observations, including the data plotted in Fi gures [2] and El suggest t hat there is a break 



at around 30 km in the KBO size distribution (e.g 


Bernstein et al. 


2004: 


Fuentes & Holman 


2008 




Fraser & Kavelaars 


2009: 


Schlichtine et al. 


2009 




Fuentes et al. 


2OI0I). This break is 



usually attri buted to collisional evolution of bodies with R < 30 km over t he age of the solar 
system (e.g. iDohnanyil Il969l : iKenyon &: Bromleyl 120041 : iPan &: Saril |2005| ). The KBO size 
distribution below radii of ~ 10 km is still poorly constrained, because KBOs of these sizes 
are too small to be detected in reflected light. They can, however, be detected indirectly, 
by stellar occultations. Recent KBO occultation surveys provide the fir st estimates for 



the abundance and upper limits of km s i zed to sub-km sized KBOs (e.g. iLiu et al.l 12008 



Schhchting et all 12009 



Bianco et al.ll2010l : ISchlichting et al.l 120121 : IZhang et al. 



20131 ). 



The work presented in this paper focuses on the size distribution of small KBOs below 
the break. We model the growth and the subsequent collisional evolution in the Kuiper 
Belt self-consistently by following the collisional evolution over 4.5 Gyrs of the whole KBO 
size distribution that resulted from runaway growth. We find that the break radius at 
i? ~ 30 km and size distribution below the break are well matched by collisional evolution 
models that use published strength laws and make testable predictions for the small KBO 
size distribution. We show that the excess mass in small planetesimals from the runaway 
growth phase leaves a permanent signature in the size distribution of small bodies that is not 
erased after 4.5 Gyrs of collisional evolution. Observations of the small KBO size distribution 
can therefore test if large KBOs grew as a result of runaway growth and constrain the initial 
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Fig. 2. — Comparison between the observed Kuiper Belt size distribution fo r objects with 
inclinations < 5°, also referred to as the cold population, as summarized in 



( 120101 ) (points), with the numerical coagulation results from lSchlichting &: Sari 



Fuentes et al. 



torn (hue) 



The error bars give the la errors on the cumulative size distribution. The observed Kuiper 
Belt size distribution above i? ~ 30 km is well matched by planet formation models of 
runaway growth. Note the break in the size distribution at i? ~ 30 km. We assumed an 
albedo of 4% and a distance of 42 AU when converting the observed magnitudes into radii. 
We note here however, that the exact choice for the value of the albedo does not affect the 
fit between the observational data and the numerical results because assuming a different 
value for the albedo would simply shift the x-axis values by a constant and this shift can be 
matched by the numerical results by letting the self-similar growth continue to larger /smaller 
sizes. 
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Fig. 3. — Comparison between the observed Kuiper Belt size distribution fo r objects with 
inclinations > 5°, also referred to as the hot population, as summarized in iFuentes et al. 



(120101 ) (points), with the numerical coagulation results from lSchlichting fc Saril (1201 ll ) (line). 



The error bars give the la errors on the cumulative size distribution. The observed Kuiper 
Belt size distribution above i? ~ 30 km is well matched by planet formation models of 
runaway growth. Note the break in the size distribution at i? ~ 30 km does not seem as 
strongly pronounced as in Figure O 
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planetesimal sizes. 

This paper is structured as follows: We describe our Kuiper Belt growth and coUisional 
evolution model in section 2. In section 3 we present our results and compare them with cur- 
rent observational constrains on small objects in the Kuiper Belt. Discussion and conclusions 
follow in section 4. 



2. Kuiper Belt Growth and Collisional Evolution Model 



2.1. Growth Model 



We use the same coagulation model as described in ISchlichting fc Saril ( 120111 ). which 
follows the mass gro wth and the coup led evolution of the veloci ty dispersion using Safron ov's 
statistical approach (jSafronovl 119721 ) . We refer the reader to ISchlichting fc Saril (120111 ) for 
the full set of equations for the growth rates of the bodies in the different mass bins and for 
the corresponding evolution of their velocity dispersions. We investigate the KBO growth 
in a single annulus centered at 40 AU from the Sun with a width of 10 AU and start the 
simulations with a total mass of about 20 Earth masses in small planetesimals. T his mass 
surface density was derived by extrapolating of the minimum mass solar nebula ( iHayashi 
198ll ) to 40 AU, after it has been enhanced by a factor of a few as required for the formatio n 



of Uranus and Neptune (e.g. iGoldreich et al.ll2004j : iDodson- Robinson fc Bodenheimerll2010l ). 



We assume that when the relative velocity exceeds the escape velocity of the larger of the 
two bodies (i.e., Vrei > VescB)^ accretion occurs and that, if the center of mass collisional 
energy of two colliding bodies exceeds the catastrophic destruction threshold, fragmentation 
takes place (see section 12^2] for details). 

In the Kuiper Belt planet formation never went all the way to completion. The growth 
was likely terminated due to the excitation of the velocity dispersion of the growing KBOs 
and small planetesimals by the formation and migration of the planets in the outer solar 
system. We model this dynamical excitation by increasing the velocity dispersion of all 
bodies in our numerical model to 1 kms~^, which corresponds roughly to the random velocity 
dispersion in the Kuiper Belt today, once objects with the size of Pluto have formed. As 
long as most of the initial mass resides in planetesimals that are about 1 km in size or larger, 
destructive collisions and fragmentation are not important until objects comparable to the 
size of Pluto have formed. This is because initially the planetesimal velocities are smaller 
than their escape velocities and even as their velocity dispersions are stirred by the growing 
KBOs, objects of at least several hundreds of kilometers in radius have to form until they can 
dynamically excite the velocity dispersion of the small planetesimals above speeds needed for 
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destructive collisionsjl] This picture changes completely once the velocity dispersions of all 
bodies are excited to 1 kms~^. From this time onward the growth is essentially terminated 
and destructive collisions lead to the onset of a coUisional cascade. We assume, although 
objects in the hot and cold population likely formed at somewhat different locations, that 
they evolve together collisionally over the age of the solar system. This assumption is likely 
valid because the same physical processes, i.e., the formation and migration of Neptune, that 
are responsible for the excitations of the KBOs' velocity dispersion are responsible for placing 
the hot population into its current l ocation. The KBO formation timescales are g enerally 



found to be less than 100 Myrs (e.g. iKenyon fc Luulll999l : ISchlichting &: Sarill201ll ). which 



suggests that the Kuiper Belt had close to 4.5 Gyrs to evolve collisionally. Observations 
of the Kuiper Belt size distribution find that the break radius and the slope of the size 
distribution below the break are the same in both the hot and cold KBO population, which 
is consist ent with the idea th at these two populations are undergoing collisional evolution 
together (IFuentes et al.l 120101 ). 



2.2. Collisional Model 

We model destructive collisions in the following way. The catastrophic destruction 
threshold, Q^, is defined as the specific energy needed to disperse the targets into a spectrum 
of individual objects such that the largest one has exactly half the mass of the original target. 
When the center of mass collisional energy of two colliding bodies, mi and m2, exceeds the 
catastrophic destruction threshold, Q}j, then the combined mass, mi + m2, is distributed 
such that one body of mass 0.5(mi + 1712) is formed and the remaining mass is distributed 
as debris over all mass bins that correspond to planetesimal sizes with m < 0.5(mi + 1712) 
according to a differential power law size distribution given by dN/ dR oc R~'^* . 

Since the Kuiper Belt consists of mostly icy bo dies with an average v elocit y dispersion 



of about 1 km s ^ , we adopt the strength law from iLeinhardt fc StewartI (|2009[ ) for ice and 



1 km s ^ impact speeds for the catastrophic destruction threshold, which is given by 

Ql = 1.3xl0'(- +0.08 ergg-i. (2) 

\ 1 cm / \ 1 cm / 

Figure m shows as a function of size and the transition from the gravity dominated regime 
{R > 0.1 km) to the material strength dominated regime (i? < 0.1 km). For comparison, 



^If initially most of the mass resided in planetesimals that were much smaller than 1 km in size, then 
the KBO growth maybe substantially different from the case investigated here, because destructive collisions 
and fragmentation will commence before bodies of a few hundred kilometers in radius have formed. 
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the catastrophic destruction threshold corresponding to the specific gravitational binding 
energy in the gravity regime and the same material strength law as before is also shown in 
Figure HI The gravitational binding energy gives an interesting absolute lower limit to the 
catastrophic destruction threshold, since bodies cannot be weaker than this. 




1 10 
Radius [km] 



1000 



Fig. 4. — Catastrophic destruct ion threshold, Qp, as a fun ction of size. The solid, blue 
line corresponds to results from iLeinhardt fc StewartI (120091 ) for ice and 1 km s~^ impact 
velocities, which corresponds to the velocity dispersion in the Kuiper Belt today. For com- 
parison, the catastrophic destruction threshold corresponding to the gravitational binding 
energy in the gravity regime and the same material strength law as before is shown as dashed 
blue line. The gravitational binding energy gives an absolute lower limit to the catastrophic 
destruction threshold, since bodies cannot be weaker than this. 



For the fragment size distribution, dN/dR oc R~'^* , we adopt q* = 3.68. This value of 
q* corresponds to the expected coUisional equilibrium size distribution, which has a power 
law index that is given by 

21 + a 



6 + a 

where a is t he exponent of R in Qn (see Equation ([2]) in the material strength dominated 
regime (e.g. 
geq = 3.68. 



Pan fc SchlichtingI |2012| ). From Equation ([3]) we find that a = —0.4 yields 



3. Results 



Combining our growth and collisional model we investigated the evolution of the KBO 
size distribution starting from various initial planetesimal sizes over 4.5 Gyrs. 



3.1. 1 km sized Planetesimals 

Figures and show the resulting KBO size distribution (solid blue line) after 4.5 Gyrs 
of growth and collisional evolution when starting from an initial planetesimal size distribution 
that consists solely of 1 km sized bodies, and from an initial planetesimal size distribution 
that has equal mass per logarithmic mass bin for bodies ranging from 0.4 km to 4 km in 
radius, respectively. For comparison, the dashed blue lines in Figures E] and [H] show the KBO 
size distribution at the end of runaway growth just before the start of destructive collisions. 
First of all it is interesting to note that the resulting small KBO size distributions do not 
follow a single power law below the break (i.e., below i? ~ 30 km) as one may naively expect. 
Instead we find that the small KBO size distribution exhibits a strong deficit of bodies 
around R ^ 10 km in size and a strong excess of bodies around 2 km in radius compared 
to abundances from a single power law size distribution spanning the range from 0.1 km to 
30 km. This deficit and excess are caused by the planetesimal size distribution left over from 
the runaway growth phase, which left most of the initial mass in small planetesimals. This 
excess mass in small planetesimals leaves a permanent signature in the size distribution of 
small bodies that is not erased after 4.5 Gyrs of collisional evolution. The resulting KBO 
size distributions shown in Figures O and M are both consistent with abundance estimates 



and upper limits from KBO occultation surveys (jSchlichting et al.ll2012l : IZhang et al.ll2013l ) 
shown in black. However, if all the mass resides solely in 1 km planetesimals initially, not 
quite enough mass is depleted in 10-30 km radius range compared to observations (see Figure 
E]). If on the other hand, we start with an initial planetesimal size distribution that has equal 
mass per logarithmic mass bin for bodies ranging from 0.4 km to 4 km in radius we find 
good agreement with the observations (see Figure [6]). 

Figure [7| shows the same collisionally evolved size distribution as in Figure M but with 
the corresponding power law indices for the different segments. KBOs with i? > 30 km 
follow a size distribution with a differential power law index g ~ 4, which is a relic from their 
formation and has not been modified by collisional evolution over 4.5 Gyrs. The power law 
index of the size distribution between 0.1 km to 30 km changes from g ~ 2 (10 km ^ R ^ 
30 km) to g ~ 5.8 (2 km < R < 10 km) and then to g ~ 2.5 (0.1 km < R < 2 km). This 
change in the slopes of the size distributions is mainly caused by the excess population of 
planetesimals that was left over from the runaway growth phase. This excess population gives 
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Fig. 5. — The small KBO size distribution after 4.5 Gyrs of collisional evolution for an 
initial planetesimal population that consisted of 1 km radii objects (blue thick line). For 
comparison, the KBO size distribution at the end of runaway growth and at the onset of 
destructive collisions is give n by the dashed blu e line. The observed KBO size distribution 
is shown by the red points (IFuentes et al.ll2010h . The black point with error bars and the 
thin black lines ranging from 0.1 km to 1 km represent the best estimate and the 95% upper 
and lower limits on the small KBO population from the HST-FGS occultation survey by 
Schlichting et al.l ( 12012| ). respectively. The thin black line ranging from 0.2 km to 20 km 



repre sents the 95% upper limit on the small KBO population from TAOS (IZhang et al. 
2018h . 



Ia 0.2 0.5 i 2 5 10"^0 50^0 
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Fig. 6. — The small KBO size distribution after 4.5 Gyrs of coUisional evolution for an 
initial planetesimal population that contained equal mass per logarithmic mass bin in bodies 
ranging from 0.4 km to 4 km in radius (blue thick line). For comparison, the KBO size 
distribution at the end of runaway growth and at the onset of destructive collisions is given 
by the dashed blue line. The observed KBO size distribution is shown by the red points 
(IFuentes et al.ll2010l ). The black point with error bars and the thin black lines ranging from 
0.1 km to 1 km represent the best estimate and the 95% upper and lower limits on th e 
small KBO population from the HST-FGS occultation survey by ISchlichting et al.l ( 12012| ). 
respectively. The thin black line ranging from .2 km to 20 km re presents the 95% upper 
limit on the small KBO population from TAOS (jZhang et al.ll2013l ). 
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rise to a very steep size distribution between ~ 2 km and ~ 10 km which grows shallower 
in time (see Figure |8]) because the excess in kilometer-sized planetesimals is being depleted 
with time. The shallow, q = 2.0, power law index between ~ 10 km and ~ 30 km is due to 
the excess population of kilometer-sized planetesimals that started to deplete the population 
of bodies between ~ 10 km to ~ 30 km. The size distribution for i? < 0.1 km takes on the 
expected equilibrium value for material strength dominated bodies as calculated in section 



2.2 fro m Equation ([3]) for the catastrophic destruction threshold from iLeinhardt fc Stewart 



(120091 ) for ice and impact velocities of 1 km s~^. The precise values of the power law index 
in the different size regimes and the exact location of the inflection points depend on the 
catastrophic destruction criterion as a function of radius and the initial planetesimal size 
distribution. For example, starting with planetesimal sizes that range from 0.4 km to 4 km 
with equal mass per logarithmic mass interval yields a smaller value for q in the 1 to 10 km 
range than starting with all the mass in 1 km-sized planetesimals (see Figures O and Figures 
[6]for comparison). However, the overall behaviour, i.e., a deficit of bodies around i? ~ 10 km 
and an excess of bodies around 2 km in radius, does not depend on the exact choices of the 
catastrophic destruction criterion and the initial planetesimal size distribution. 

Figure [H] displays the time evolution of the small KBO size distribution. The differential 
power law indices between ~ 10 km and ~ 30 km and between ~ 2 km and ~ 10 km 
become shallower with time. The decrease in the power law index between ~ 2 km and 
~ 10 km is due to the fact that the excess population of planetesimals that was left over 
from the runaway growth, which gave rise to a very steep size distribution between ~ 2 km 
and ~ 10 km (dashed blue line in Figure E]), is being depleted by destructive collisions with 
time. The decrease in the power law index between ~ 10 km and ~ 30 km is due to the 
excess population of kilometer-sized planetesimals that start to deplete the population of 
bodies between ~ 10 km t o ~ 30 km. I.e., the catastrophic destruction threshold from 



Leinhardt fc StewartI (120091 ) yields, for a velocity dispersion of 1 km s ^, that 10 km sized 



bodies are typically destroyed by bodies ~ 1 km in radius and 30 km sized bodies are 
typically destroyed by bodies that are ~ 8 km in radius. The power-law index below ~ 2 km 
evolves to g ~ 2.5 and remains close to constant from then onwards. 



3.2. 10 km sized Planetesimals 

Figure M shows the resulting small KBO size distribution after 4.5 Gyrs of collisional 
evolution for an initial planetesimal population that consisted of 10 km radii bodies (solid 
blue line). The small KBO size distribution is inconsistent with the observed size distribution 
of KBOs with radii ranging from 10 km to 100 km (red points) and with upper limits from the 
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Fig. 7. — The same small KBO size distribution after 4.5 Gyrs of collisional evolution as 
shown in Figure [6] but plotted with the corresponding differential power law indices for the 
different segments of the size distribution. The deficit around 10 km results from an excess 
of ~ 1 km planetesimals at the onset of the collisional evolution. The size distribution for 
i? < 0.1 km takes on the expected equilibrium value for material strength dominated bodies 
as calculated in section 2.2. The size distribution above i? ~ 30 km remains unchanged by 
collisional evolution over the age of the solar system and is therefore primordial. 
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Fig. 8. — The time evolution of the small KBO size distribution for an initial planetesimal 
population that contained equal mass per logarithmic mass bin in bodies ranging from 0.4 km 
to 4 km in radius. The thin blue line corresponds to the KBO size distribution at the end of 
runaway growth before the onset of destructive collisions; the yellow, dotted line corresponds 
to 100 Myrs; the green dot-dashed line to 1 Gyr; the solid, blue line to 4.5 Gyrs, and the 
dashed purple line to 10 G yrs of collisional evo lution. The observed KBO size distribution 
is shown by the red points ( iFuentes et al.ll2010l ). The thin black lines ranging from 0.1 km 
to 1 km represent the 95% upper and lower limits on the small KBO population from the 
HST-FGS occultation survey by ISchlichting et al.l ( 120121 ) and the thin black line ranging 
from 0. 2 km to 20 km rep resents the 95% upper limit on the small KBO population from 
TAOS jzhang et al.lboisL 
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TAOS KBO occultation survey ( iZhang et al.ll2013l ). Even if we assume that large KBOs are 
only held together by their own gravity (blue dotted line in Figure [9]), which is an absolute 
lower limit on their strength, because bodies can't be weaker than this, we find that we 
cannot match the observed KBO size distribution. We also started with initial planetesimal 
populations that contained equal mass per logarithmic mass interval between 1 km and 
10 km in radius and were still unable to find a reasonable agreement between the resulting 
small KBO size distribution and the observations. This leads us to conclude that the Kuiper 
Belt did not form via coagulation from an initial planetesimal population that contained 
most of the initial mass in planetesimals that were 10 km in radius or larger, because not 
enough of the initial planetesimals can be ground down over the age of the solar system to 
match observations. 



4. DISCUSSION AND CONCLUSIONS 

We studied the size distribution of small KBOs by modeling self-consistently the growth 
and the subsequent coUisional evolution over 4.5 Gyrs in the Kuiper Belt and arrive at the 
following results: 

• The Kuiper Belt size distributions of the cold and hot population for radii > 30 km are 
primordial and can both be well fit by the resulting size distributions from planet formation 
models of runway growth (see Figure [2] and [3]) with the notable difference that the largest 
bodies in the hot population grew to larger radii than in the cold population, which is 
consistent with the idea that the hot population formed at smaller semi-major axis compared 
to the cold population. 

• The break radius at i? ~ 30 km and size distribution below the break are well matched 



by collisional evolution models that use published strength laws ( iLeinhardt fc Stewartll2009l ) 
and start with resulting size distributions from runway growth and follow the collisional 
evolution in the Kuiper Belt over 4.5 Gyrs. 

• Compared to a single power law size distribution that would span the whole range 
from 0.01 km to 30 km, we find in general a strong deficit of bodies around i? ~ 10 km and 
a strong excess of bodies around 2 km in radius. This deficit and excess are caused by the 
planetesimal size distribution left over from the runaway growth phase, which leaves most of 
the initial mass in small bodies. This excess mass in small planetesimals leaves a permanent 
signature in the size distribution of small bodies that is not erased after 4.5 Gyrs of collisional 
evolution. Future KBO occultation surveys, which probe the small KBO size distribution, 
can therefore test if large KBOs grew as a result of runaway growth and constrain the initial 
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Fig. 9. — The small KBO size distribution after 4.5 Gyrs of collisional evolution for an initial 
planetesimal population that consisted of 10 km radii bodies (blue thick line). The resulting 
size distribution assuming that large KBOs are solely held together by their own gravity is 
shown by the dotted blue line. For comparison, the KBO size distribution at the end of 
runaway growth and at the onset of destructive collisions is given by the dashed blue line. 
The observed KBO size distribution is shown by the red points. The black point with error 
bars and the thin black lines ranging from 0.1 km to 1 km represent the best estimate and 
the 95% u pper and lower limits oii the small KBO population from the HST-FGS occultation 
survey by 



Schlichting et al.l (|2012[ ). respectively. The thin black line ranging fr om 0.2 km to 



20 km represent the 95% upper limit on the small KBO population from TAOS (jZhang et al. 



20131 ). The resulting KBO size distributions are inconsistent with the observed KBO size 
distribution for bodies with radii ran ging from 10 km to 100 km and with upper limits from 
the TAOS KBO occultation survey (jZhang et al.ll2013l ). This results holds true even if we 
assume that large KBOs are only held together by their own gravity. We conclude that 
the Kuiper Belt did not form from 10-km sized planetesimals by coagulation, because not 
enough of the initial planetesimals can be ground down over the age of the solar system to 
match observations. 
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planetesimal sizes. 



The observed KB O size distribution derived by iFuentes et al.l ( 20101) by cqmbin- 



2008 



ing various KBO surveys (Chianjg fc Brown 19991: Gladman et al 



Allen et al.lbooilBernstein et al.ll2004j: 



Fuentes et al. 



2009 



Petit et al.ll2006 



Eraser fc Kavelaars 



cent optical KBO occupation surveys ( ISchlichting et al. 



2009 



2001 



Eraser et al. 



Fuentes et al. 



2012 



2008 



ruiillo et al. 



2001 



Fuentes fc Holman 



2010 ) and results from re- 



Zhang et al.l l2013l ) are best 



matched by an initial planetesimal population that contained about equal mass in bodies 
ranging from 0.4 km to 4 km in radius. In addition, the resulting KBO size distribution after 
4.5 Gyrs of coUisional evolution is also consistent with upper limits from KBO occultation 
surveys at x-ray wavelengths that probe o bjects ranging from ~ 30 m to 300 m in radius 
J Jones et al.ll2008[ Ichang et allboill . l2012h . 



• The observed KBO size distribution for R > 10 km cannot be matched if most of the 
initial planetesimal mass resided in bodies that were 10 km in radius or larger, because their 
resulting size distribution cannot be sufficiently depleted over 4.5 Gyrs to match observations. 
We conclude from this that the Kuiper Belt did not form by coagulation from an initial 
planetesimal size distribution that contained most of its mass in planetesimals with radii of 
10 km or larger. 

There are several further interesting things to note here: 

Since the excess mass in small planetesimals from the runaway growth phase leaves 
a permanent signature in the size distribution of small bodies that is not erased after 4.5 
Gyrs of coUisional evolution, future KBO occultation surveys will be able to test whether 
large KBOs grew as a result of runaway growth from an initial planetesimal population 
consisting of bodies ranging from a few hundred meters to a few kilometers in size. The 
small KBO size distribution therefore offers the opportunity to observationally constrain the 
initial planetesimal sizes from which planets form, which re mains one of the major open 
questions in planet formation theory (jChiang &: Youdinll2010l ). 



The resulting small KBO size distributions that we find a ll contain enough bodie s to 
satisfy the required supply rate for the Jupiter family comets (IVolk &: Malhotral 120081 ) . If 
the Kuiper belt formed by coagulation from km-sized planetesimals then there should be an 
excess of about a factor of 1000 of small comets with initial radii (i.e., before any mass loss 
or break up occurs) of 2 km compared to 10 km. 

Because the comet size distribution has likely been modified by mass loss and break up 
of the cometary nuclei, the size distribution of centaurs should provide a more reliable probe 
of the KBO size distribution between 1 to 10 km in radius. Unfortunately, although about 
two hundred centaurs are currently known with sizes ranging from about 100 km to 1 km. 
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no well characterized survey has been carried out to date that would allow the derivation of 
a de-biased centaur size distribution and therefore probe the small KBO size distribution. 

Finally, it is very intriguing that there is a striking similarity between the small KBO 
size distribution that we find after 4.5 Gys of coUisional evolution and the reconstructed 



impa ct size distributi o n from the cratering records on the Saturnian satellites ( iMinton et al. 



20121). iMinton et al.l ( 120121 ) find that the cratering size distribution of the old terrains of 
Dione, Hyperion, lapetus, Mimas, Phobe, Rhea and Tethys can be explained by a single 
impactor population that follows a size distribution with differential power laws indices of 
g = 4 for i? > 30 km, q = 2.0 for 10 km < R < 30 km, q = 4.2 for 1 km < R < 10 km, 
q = 2.6 for 0.1 km < R < 1 km and q = 3.7 for 0.01 km < R < 0.1 km. These values are 
in remarkably good agreement with the power laws indices that we find for the small KBO 
size distribution and that are shown in Figures [71 The only notable di fference between our 
small KBO size distribution and the results from IMinton et al.l (120121 ) seems to be in the 
range from ~ 1 km to ~ 10 km for which we find a steeper size distribution with power law 
index q = 5.8. The similarities between our results for the small KBO size distribution and 
the reconstructed impactor size distribution suggest that the impactors that bombarded the 
Saturnian moons originated from the Kuiper Belt. 



We thank David Jewitt for his comments and suggestions. 
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